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Abstract 



In this paper, we are concerned with the controllability of a chemotaxis system of parabolic- 
i—i .... 
^. , elliptic type. By linearizing the nonlinear system into two separated linear equations to bypass 

the obstacle caused by the nonlinear drift term, we establish the local null controllability of the 

p — ■ original nonlinear system. The approach is different from the usual way of treating the coupled 

•/"") ■ parabolic systems. 

Keywords: Local null controllability, chemotaxis system, parabolic-elliptic type, Kaku- 
tani's fixed point theorem. 

AMS subject classifications: 93B05, 93C20, 35B37. 



_<3 ■ 1 Introduction and main result 

In this paper, we are concerned with a controlled initial-boundary value system of parabolic-elliptic 
type 

8 t u = V ■ (Vu - xuVv) + l w / in x (0, T), 

= Av - -yv + 5u infix (0, T), 

(1.1) 

d y u = 0, d y v = on 89 x (0, T), 

u(x,0) = Uq(x) X £ 9, 
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where and henceforth u and v are shorthands for states u(x,t) and v(x,t) at time spacial position 
iEf! and time t > 0, dt = d/dt, d v = d/du is the directional derivative along the outward unit 
normal vector v on d£l, l w stands for the characteristic function of w, / is the control function, no 
is the initial value, %, 7 and 5 are given positive constants. 

The system (|l.ip without the control (i.e., / = 0) is a simple chemotaxis system, which was 
addressed by Keller and Segel p3] as a model to describe the aggregation process in slime mold 
morphogenosis, assuming that the cells emit directly the chemoattractant which is immediately 
diffused. The unknown function u then denotes the cell density, whereas v represents the concen- 
tration of the chemoattractant. The validity of Keller-Segel's chemotaxis system is supported by 
some experiments on the Escherichia Coli bacteria and other interesting physical interpretations. 
In fact, it has been extensively involved in many medical and biological applications as well as some 
relevant areas such as ecology and environment sciences. Moreover, because the model has a rich 
structure from mathematical point of view, it deserves to be challenged. Actually, some special 
but interesting cases have been studied such as the aggregation, the blow-up of solutions, and the 
chemotactic collapse. Some significant results have been achieved from different perspectives. We 
refer to [TT] (also |12j ) for a survey where a quite complete bibliography on the topic is included. 

In this paper, we study the Keller-Segel system from the controllability point of view. We say 
that the system (jl.lh is locally null controllable at time T, if there exists a neighborhood of the 
origin such that for any initial data uq belonging to this neighborhood, the solution (u, v) of (|1.1|) 
produced by corresponding control function / satisfies u(x,T) = for almost all x G CI, where the 
neighborhood and the control function space will be specified later. Here, we consider the local 
null controllability instead of the exact null controllability. The reason being interested is that the 
solutions of Keller-Segel system may blow up in either finite time or infinite time, which is shown 
in [10] for the 3-d case. The 2-d case is even more interesting and attractive. Actually, it has been 
found for the 2-d case that the solution exists globally in finite time when the mass of the initial 
value is less then a threshold value, while the solution will blow up either in finite or in infinite 
time when the mass of initial value is larger than the threshold value (see, e.g., |13|). 



The study of the controllability for parabolic equations has been thriving in the past decade, see 
for instance [51 [B] and the references therein. Among them, a special interest is on the controllability 
of coupled parabolic systems. For coupled systems, the most practical situation is to impose the 
control force on one equation, which has attracted intensive attention in the last few years. We 
refer to a survey paper [2] wherein abundant references are provided. 

However, to the best of our knowledge, very few results are available to the control problems 
of the Keller-Segel system (jl.ip where a parabolic equation is coupled with an elliptic equation 
through a drift term. Very recently, the null controllability of a kind of nonlinear parabolic-elliptic 
system of the following is considered in [7j : 

d t y-Ay = F(y,z) + l UJ f in(lx(0,T), 
-Az = f(y,z) inOx(0,T), 



where F(y,z) and f(y,z) are nonlinear terms. In [16], an optimal control problem of the system 
(jl.ip with the control to be imposed on the second equation is considered. Our previous work 
[9] is the first work that considers the local exact controllability of a type of Keller-Segel system 
where a parabolic equation couples with another parabolic equation. In the present work, we 
attempt the controllability of system (II. ip . which is probably the first work for this system. In 
the system (jl.ip . since the drift term — %V • (uVv) destroys some good properties of the diffusion 
operator which ensure the regularity of the system, much more mathematical difficulties than the 
aforementioned coupled parabolic systems are caused. These include the regularity, the estimation 
of the "observability inequality", and many others. The usual way to establish the controllability of 
a nonlinear system is to linearize the nonlinear system into some coupled linear ones. Then, by the 
combination of the controllability result of the linearized system and some fixed point results, one 
is able to establish the controllability of the nonlinear system. We refer to the typical works E] 
for the approach of this kind. In this paper, we investigate, however, the controllability of system 
(|l.ip in a different way motivated intuitively by the special mathematical structure of system (jl.ip . 
We may decompose this nonlinear system into two separated linear systems: one is a controlled 
parabolic system, and another one is an irrelevant elliptic equation. In such a way, we can bypass 
the obstacle caused by the nonlinear drift term. This technique would be useful for other coupled 
systems like drift-diffusion equations from the semiconductor device. 

Throughout the paper, C Wl N (N > 1) stands for a bounded domain with smooth boundary 
dfl, to is a nonempty open subset of 0, and T > 0. Q = Q x (0, T), £ = dQ x (0, T), and 
Qu = u) X (0, T). The norms of the usual Lebesgue function spaces L P (Q) and L P (Q) are denoted 
by | • |p and || • || p , respectively. W s > q (n), Wq' X (Q), C 2 > l {Q) and C a {tt) (s, a > 0, 1 < q < oo) 
represent the usual Sobolev spaces (see, e.g., [IS]). When q = 2, H m (£l) = W m,2 (Q), m G N. In 
addition, 

W(0,T) = {y\y G L 2 (0, T; H 1 ^)), 8 t y G L 2 (0, T; H 1 ^)*)} 

is equipped with the graph norm \\y\\ W (o,T) = IMIz^CojVff 1 ^)) + \\dty\\L 2 (o,T;m(n)*), where H 1 ^)* 
denotes the dual space of H 1 (Q) and their duality product is denoted by (•, •}. We also use C to 
denote a positive constant independent of T, which have different values in different contexts. 

Definition 1. A pair of functions (u,v), with u G W(0,T) n L°°(Q) and v G L 2 (0,T; H 1 ^)) n 
L°°(Q), is said to be a weak solution of (jl.ip if for every ip G L 2 (0, T; H l (Q)), the following 
identities hold true: 

/ (d t u,(p)dt+ // [(Vu - xu^v) ■ Vip - luftp] dxdt = 0, 
Jo J Jq 

/ / [Vv ■ Vtp + (jv - 5u)ip] dxdt = 0. 

Now, we are in a position to state the main result of this paper. 
Theorem 1.1. Let T > 0. For any initial value uq satisfying 

< e - cl (l+T+±) ; (L2) 



where c\ > is a constant independent of T , there exists a control f G L°°(Q U] ) such that system 
(jl.ip admits a solution (u, v) satisfying 

u G C([0, T]; L 2 (n)) n L°°(Q), u G L 2 (0, T; fl 4 ^)) D £°°(Q), 

and u(x, T) = /or almost all xGfl. 

Remark 1.1. Theorem 11.11 gives an explicit representation of the initial data with respect to time 
T. It particularly shows that the shorter of the terminal time T is, the smaller for the initial value 
for the null controllability of system (|1.1|) . In addition, under the assumption of Theorem 11.11 
limsup t _+ r - ||f(-,t)||2 = 0. 

We proceed as follows. In section 2, we give some preliminary results. Section 3 is devoted to 
the null controllability of a scalar parabolic equation, for which the L°°-control is obtained and its 
estimates with respect to time T are also established. The proof of Theorem 11.11 is presented in 
section 4. 

2 Some results for linear equations 

In the sequel of the paper, we need some regularity results for linear equations for both parabolic 
and elliptic types. We first consider the well-posedness of the linear elliptic equation followed by 

= Av — ryv + Sri in Q, 

(2.1) 

d u v = on dQ, 

where 7 and 5 are positive constants. The result of Proposition 12 . 1 1 is brought from pp. 
Proposition 2.1. For any rj £ L p (Q),p > 1, Eq. (12. ip admits a unique solution v £ W 2,P (Q) with 

\\v\\ W 2, P(Q) <C\rj\ p . 
Next, we consider the parabolic equation 

d t u = Ait - V ■ (Bu) + F in Q, 

d v u = on E, (2.2) 

u(x, 0) = uq(x) x G 0,. 

Proposition 2.2. Let B G L°°(Q) N with B ■ v = on E, F G L°°(Q), and u G L°°(ft). TTten 
Eq. (|2.2p admits a weak solution u G L°°(Q) with 

||«||oo<e c « (|u |oo + ||i r ||oo), (2-3) 

where C = C(Q) is a positive constant depending on 0, and go is given by 

Q0 = (l + \\B\t)(l+T). (2.4) 



A similar inequality (|2.3p could be found in [15], but here we improve the estimate so that it 
depends on time T explicitly. To this end, we need the following lemma (see |15|. Lemma 5.6, p. 

95]). 

Lemma 2.1. Suppose that a sequence Y s , s = 0, 1,2, • • ■ of nonnegative numbers satisfy a recursion 
relation 

Y s+1 <cb s Y s 1+e , S = 0,l,2,--- 

with some positive constants c,e and b > 1. Then Y s — > as s — >■ oo provided that 

Y < C~H~e*. 

Proof of Proposition 12.21 Let (u — k) + = max{u — k, 0} and Ak(t) = measjx € U\u(x, t) > k} 

for t G [0,T], where 

k>K = \\F\\ oo + \u \ oo . (2.5) 

Multiplying by (u — k) + the both sides of (|2.2p . we get, by integration by parts and the Holder 
inequality, that 

— / \(u — k)+\ dx + / |V(u — fe)+| dx 
dt Jn Jn 

< ||S||^ / u 2 cfc 



/|(u-fc) 

<(2|| J B|| 2 0O + l)( f \( u -k) + \ 2 dx+ f 
\Jn JA k 



A k (t) Jn JA k (t) 

k 2 dx 



At) 
From Gronwall's inequality, it follows that 

\{u-k) + \ 2 dx + f [ \V(u-k) + \ 2 dxdt<e Cgo [ [ k 2 dxdt (2.6) 

n Jo Jn Jo J A k (t) 

for all t € [0, T], where and in what follows qq is given by (|2.4p . On the other hand, by Proposition 
1.3.2 of [6], 

inky < c(n)(i + t)to [|«||^ (Q) (2.7) 

for any v € V^Q). Here V 2 {Q) = L°°(0,T; L 2 (0)) n L 2 (0,T; H l (tt)) is endowed with its graph 
norm. Then, (|2.6p together with (|2.7p gives 



||(tt-fc)+ 11 2(^+2) <e° eo / / rdxdi. (2.8) 

" Jo JA k (t) 

Let <p(k) =meas{(x,t) G Q|n(x,t) > A;}. Then for any h > k, we get, from (|2.8p . that 

(/i - k) 2 tp(h)^ < ||( n _ fc) + || 2 2 ii v ± 2) < e c e°<p(k)k 2 , 

N 

which then gives 

2(iV+2) 

^ h ) ^ e ° go (jr^k) v(k)^- (2-9) 



Next, set Y s = ip(k s ), k s = M(2 — ^), and put h = k s +i and k = k s in (j2.9[) to get 

Y s+1 <c4 T (2 T ) s Y s 1+e , 

where r = 2(N + 2)/iV, e = 2/JV, and c = e Cgo . By Lemma EH we get <p(2M) = provided that 

Y = <p(ko) = <p(M) < (cA T )—c (2 r y^, (2.10) 

for some positive real number M. To determine the value of M, let m > 1 be an integer and 
M = itlKq, where Kq is given by (|2.5[) . Put h = M = uiKq and k = Kq in (|2.9[) to get 

y(M) < c ( — !— -^ ^(i^o) 1+e < c ( — ^-^ T 1+e (measft) 1+e . (2.11) 

\m — 1/ V m — 1/ 

Now, to get </j(2M) = 0, we need to choose a proper m such that (|2.1U|) holds. Combining (|2.10p 
and (|2.1ip . we only need the integer m to be such that 

l+g l+g ,-»!+£ -+ 1 

m > 1 + c st T t (measi2) t 2 e 7 1 . 
Hence ip(2M) = (p(2mKo) = gives 

u<2mK Q <e c ^{\\F\\ oo + \u \ oo ). 



In a similar argument, we can also get the other half part of (|2.3p for —u. This completes the 
proof. □ 

3 Null controllability of a linear parabolic equation 

In this section, we consider the null controllability of the linear parabolic equation 

d t u = An - V • (Bu) + l u f in Q, 

d v u = on E, (3.1) 

u(x,0) = uq(x) x G £1. 

Theorem 3.1. Let T > 0, and B G L°°(Q) N with B ■ u = on E. For any u G L 2 (0), £/iere 
exists a control f G L°°{Q iJ j) such that the solution u of system fj3. 1 j) corresponding to f satisfies 
u G W(0, T) and n(x, T) = /or ie!l almost everywhere. Moreover, the control f satisfies 

HWIL^kl^ (3.2) 

where 

« = (l + ||B||^)(l + r) + i (3.3) 

To prove Theorem l3.lt we need to establish a type of "observability inequality" for the following 
adjoint equation of f)3. lj) : 

f -d t (j) = A(j> + B- V(j> inQ, 
d v cp = on E, (3.4) 

(/>(x,T) = c/) T (x) xeQ, 

6 



where (j) T G L 2 (f2). By [SJ Lemma 1.1], there is a function /3 G C 2 (Q) such that /3(x) > for all 
ieU and /3|an = 0, 1 V/3(x)j > for all x G \ w. For A > 0, set 

e A/3 e A/3 _ e 2A||/3|| c{TT) 

^ = ^r37)' a= t(T-t) ■ (3 ' 5) 

We then have a Carleman inequality stated in Lemma 13.11 (see [8] ) . 

Lemma 3.1. There exists a constant Xq = Xq{Q,u) > 1 such that for all A > Ao and s > 

7(A)(T + T 2 ), 

/"/" [s^|Vy| 2 + {sip) 3 \y\ 2 ] e 2sa dxdt 
< C [I e 2sa \d t y ± Ay\ 2 dxdt + If (s^) 3 e 2sa \y\ 2 dxdt (3.6) 

JJq J Jqu, 

for all y G X = {£ € C 2,1 ((5)|9i/C = on £}, where 7(A) is given by 

7 (A) = e 2A|l/3|| c(«). (3.7) 

Proposition 13.11 is an "observability inequality" for the adjoint equation (13.41) . 

Proposition 3.1. Let 5q G (1,2). ITien there exist positive constants A and s such that for all 
T > 0, 4> T G L 2 (0), i/je solution <f) of system (|3.4p satisfies 



where K is given by ([3. 3D . 



■,0)\l<e CK ll e 5osa \0\ 2 dxdt, (3.8) 



Proof First, by Lemma |3.1[ there exists a positive constant Ai = Ci(il,o;)(l + K-Bjl^) satisfying 
7(Ai) > Ai > 1 such that for any A > Ai, s > 7(A) (T + T 2 ) and (p T G L 2 (£l), the associated solution 
(j) to (|3.4j) satisfies 

// (s<p) 3 \(j)\ 2 e 2sa dxdt <C II (sip) 3 |0| 2 e 2sa dxdt, (3.9) 

where C = C(fi,w), and 7(Ai) and 7(A) are given by (|3.7p . By 

j t (>" 2 ~ f m 2 .) > 0. 
This gives, for any t G (0, T], that 

|0(-,O)|^< e "^|<X,i)|^ (3.10) 

for all t G (0,T]. Integrate both sides of (EHUD over [T/4, 3T/4] to give 



l-BII 2 T 1 // .\|2 



: J ,T 



i o)| 2 <l e ii^iLT/ 4 



4 



Since (s<p)~ 3 e~ 2sa < e Cs / T2 in VL x [T/4,3T/4], inequality (|22D then follows from ([MD and (IXTTT) . 
with A and s taken as A = C(l + \\B\^) > Ai and s = C(l + ||-B||^)(T + T 2 ). □ 



Proof of Theorem 13.11 Let e > 0. We set 

J £ (u,f) = \tt \f\ 2 e- 5osa dxdt + y t \u(x,T)\ 2 dx 



'Q 

and consider the extremal problem 



inf J £ {u,f), 

(u,f)GU 



where U is the totality of (u, f) € W(0, T) x L 2 (Q) solving (|3.ip . The existence of an optimal 
pair (f £ ,u £ ) to the above extremal problem follows from a standard argument. By the maximum 
principle (see 0]), we get the optimality system for this problem as follows: 

-d t (f) e = A(b e + B- Vcj) £ in Q, 
d v <j> e = on E, (3.12) 

<j) e (x,T) = -\u e {x,T) xen- 

d t u e = Au £ - V • (Bu £ ) + l w f £ in Q, 

d v u e = on S, (3.13) 

u £ (x, 0) = uq(x) x G £1; 
f £ - l^e 50 ^ = 0. (3.14) 

Note that we can choose s and A such that the observability inequality (13.8|) and 



W (A) = e - A ^Hc(n) < S -1 (3.15) 

hold. Furthermore, by (|3.12p . (|3.13p . (|3.14p . and Proposition 13.11 we get 

ft \0 £ \ 2 e Sosa dxdt+- f \u £ (x,T)\ 2 dx<e c ' K \u \ 2 2 . (3.16) 

Eq. (|3.14p and (|3.16p then lead to ||l w / £ || 2 < e CK \uq\ 2 , which means that the controls f £ can be 
taken in L 2 space. 

Now, we show that f £ can actually be taken in L°° space. To this purpose, we apply a so-called 
bootstrap method in |17| (see also [5]). Firstly, set ao = mm-^a. Then the following inequalities 
can be easily verified: 

«o < a < rr^Tu < °' ( 3 - 17 ) 

1 -f- 0J\A) 

where w(A) is defined by (I3.15J) . Secondly, let r be a sufficiently small positive constant and let 
{tj}- be a finite increasing sequence such that < Tj < r, j = 0, 1, . . . ,M,tm = t. Let fe}j =0 
be another finite increasing sequence such that po = 2,pm = oo, and 

AT \ (\ 1 \ 1 . 

- + 1 +1>-, i = 0,l,...,M-l. 3.18 

2 ) \pi p i+ ij 2 

For each i, i = 0, 1, . . . , M, define 

z 4 (x,t) = e( s+ ^) Q0 ^(x,r-t), 
Fi(x,t) = [a t (e( s+ ^ Q0 )]^(x,T-i), 

B(x,<) = B(x,T-t). 



Then, the adjoint equation (|3.12p is transformed into the following initial-boundary problem for 
every i = 0, 1, . . . , M: 

d t Zi - Azi = B -Vzi+Fi in Q, 

d vZi = on S, (3.19) 

Zi(x,o) = ie!l. 

Let {S(t)}t>o be the semigroup generated by the Laplace operator with Neumann boundary con- 
dition. It follows that (see, e.g., [3]) 



\S(t)uL < Cm{t) 2 Vp V \u\ 



(3.20) 



for all u £ L p (il), t > 0, and 1 < p < q < oo, where m(t) = min{l, t}. Note that the solution Zi of 
(|3.19|) can be represented as 



Zi (-,t) = / S(t-s)[B-Vzi + Fi)(-,s)ds, » = 0,1,... ,M. 



(3.21) 



Applying the estimates ()3.20p to (j3.21|) . we have 



\zi{;t)\ V4 <C J m(t-s) 



2 \Pi-i Pi 



B-Vzi + Fi)(-,s) ds, i = 0,l,...,M. (3.22) 
Pi-i 



With (13.18J) . we apply Young's convolution inequality (see, e.g. [3]) to the right-hand side of (|3.22ft 
to get 

llr-ll < f> c '( 1+T ) Ml Rll llVr-ll J- II F. II | ('3 93'! 

\\Zt\\ p .^e [\\ J:S \\ 00 \\vz l \\ p ._ l -t\\f l \\ p ._ 1 j . [6.Z6) 

On the other hand, by a standard energy estimate applied to (I3.19p . we can get the following 
jj>i-i -estimate for Z{\ 



+ l|Vzi|L <e CK \\Fi\ 

' 1 1 

By the definition of Fi and ()3.17|) , we have 



m\ P ^< CT \\^K-, 



(3.24) 



(3.25) 



Combining ([333]) . (ET24"|) . and (pT25l) then, we obtain 

llz-ll < e° K \\z- ill 
||^||p. ^ e ||z l _i|| p ._ i , 

This iteration inequality from to M produces 

\\zm\\ Pm <e CK ||z || 2 - 



(3.26) 



Since pm = oo, it follow from the definition of Zq, ()3.16|) . and ()3.26|) that H^mIIoq < e° K ||ito|| 2 ! that 
is, ||^ £ e( s+T ) ao || < e CK \\u \\ 2 . By (gUD , we get 



e [-s(5o-l-w(A))+T(l+ £ j(A))]a 1 ^j £ 



<e CK ||n || 2 , 



(3.27) 



where ui(X) is given by (|3,15p . By choosing r small enough so that — s (Sq — 1 — w(A))+r(l+u;(A)) < 
0, we conclude from (|3.27[) that 

HWelloo^e^Huolla- (3-28) 

This shows that the controls f e can be taken in L°° space. 

Finally, by (|3.28p . we can extract a subsequences of {/ £ } £ >o, still denoted by itself, such that 
lu/e ~~ > lw/ ^ L°°{Q) weakly in L 2 (Q) as e — ► 0. Denote by u £ the solution to the system (|3.13p 
associated to f £ . By virtue of Proposition 12.21 { u e}e>o is uniformly bounded in W(0,T). Thus, we 
can extract a subsequence of {u £ } £ >o, still denoted by itself, such that u £ — > u weakly in W(0,T) 
for u G W(0,T) C C([0, T]; L 2 (Sl)). Such a it is the weak solution of (|3.ip corresponding to /. In 
addition, by (|3.16p . u(x,T) = almost everywhere in SI. This completes the proof. □ 

4 Proof of Theorem 11.11 

Let K = {£ G L°°(Q)| H^l^ < 1} D L°°(0, T; If (SI)) C £ 2 (Q), p > max{iV,2}. For every £ G if, 
consider the following two linear equations: 

= Av(;t)-Tv(;t) + 6£(;t) in SI, 
d u v(-,t) = on dfi 

for almost every t G [0, T] and 

<9 t u = An - V • (Bu) + l u f in Q, 
d u u = on £, 

u(x,0) = Uo(x) x G 0, 

where .B = xVu ■ I n what follows, we denote by shorthand v* = v^(-,t) the unique solution 
of equation (14. ip corresponding to £(•,£) for £ G [0,T]. First, by Proposition 12.11 we see that 
v € £ L°°(0,T;W 2 ' P (SI)) for p > maxjiV, 2} provided that £ G K . Hence, the embedding theory 
between Sobolev spaces for p > N (see, e.g., [6]) then implies 

B = X W G L°°(0, T; W^iSl))* C L°°(0, T; (7(0))* with B ■ v = on S, 

and in addition, 



Moo A. 



Vw ? 



< Clieiloo < C. (4.3) 

L°°{0,T;W^P{n)) 



Thus, we can define a linear continuous operator <£ from K to L°°(0, T; (7(0))^ C L°°(Q) N by 

$(f) = B = X V^,V^GK 

Second, by Theorem 13. 1| we see that for each B G L ao (Q) with B'!/ = 0onE, there exists a pair 
(u,f) G £ 2 (Q) x L°°(Q) that solves system ([42]) with u(x,T) = for almost all x eSl. Moreover, 
the control / satisfies §3?2$. By (|Q1) . 

HWIL^e^KIs, (4.4) 

10 



where 

Ko = l + T + i (4.5) 

and in the sequel, C is a positive constant independent of time T. By (|2.5p in Proposition 12.21 and 
(|4.4[) . we have the following estimate: 

\H\w(0,T) + Hloo < eCK ° Moo ■ ( 46 ) 

We then define a multi- valued mapping \& : L°°(Q) N — >• 2 L (^ by 



*(B) = u£ L 2 (Q) 



3/ G L°°(Qu) satisfying (|4.4p such that it is the solution 
of (J4.2J) corresponding to / and B, and u(x, T) = a.e. in f2 



where 2 L (Q> stands for all subsets of L 2 (Q). Since both operators $ and Vl/ are well defined, which 
is guaranteed by Proposition 12.11 and Theorem 13. 1\ we let 

A = fo$:i(cL 2 (Q)4 2 t2W) . (4.7) 

Now, we apply Kakutani's fixed point theorem (see [H p. 7]) to the map A to prove Theorem 1 1.11 
Indeed, it is clear that K is a convex subset of L 2 (Q). By Proposition ^. II and Theorem l3. II again, for 
any £ G K , A(£) is nonempty and it is also convex due to the linearity of the equations. Moreover, 
from (|4.6p . it follows that for each £ G K, A(£) is bounded in W(0, T), and hence a compact subset 
of L 2 (Q) according to the Aubin-Lions lemma (see [H p. 17]). 

We claim that A is upper semi-continuous. Indeed, let {Cn}^Li be a sequence of functions in K 
such that 

in — >■ £ strongly in L 2 (Q) as n — ► oo. (4.8) 

For every n, let B n = 3>(£ n ) = xVi> n and take u n G A(£ n ) = \l/(.B n ), where i> n solves 



(4.9) 



f 



= Av n (-,t) - jv n (-,t) + 5U(-,t) in n, 
d u v n (-,t) =0 on d£l 

for almost every £ G [0, T] and u n solves 

d t u n = Au n - V • {B n u n ) + l w /„ in Q, 

d u u n = on S, (4.10) 

u ra (x,0) = no(a;) i£(!, 

with u n (x,T) = for almost all xell. Moreover, the control f n satisfies 

HWnlloo^e^Ma- (4.H) 

To show that A is upper semi-continuous, it suffices to prove that there exist a subsequence of 
{«„}"=! such that it converges strongly to an element of A(r/) in L 2 (Q) topology. 

In what follows, we do not distinguish the sequence and its subsequence by abuse of notation. 
First, the estimate (|4.1ip enables us to obtain a function / G L°°(Q) and a subsequence of {/ n }^Li 5 
such that 

lw/n -> lw/ weakly in L 2 (Q); weakly* in L°°(Q) as n -> oo. ( 4 -12) 
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By (|4,lip and Proposition 12.21 u n satisfies (|4,6p ; that is 



*n|IW(0,T) 



<e CKa \u \ 2 . (4.13) 



Applying the Aubin-Lions lemma again, we get a u G W(0,T) n L°°{Q) and a subsequence of 
{^n}^Li such that 

■u n — > u weakly in W(0, T); strongly in L (Q), as n — > oo. (4-14) 

Furthermore, by the strong convergence of {u n }^ =1 in L 2 (Q), we can extract a subsequence of 
{un}^Li ( see [13 Lemma 2.1, p. 72]) such that 

u n — > u almost everywhere in Q as n —> oo. (4-15) 

On the other hand, by Proposition 12. 1\ for each n and p > 1, it holds that 

IM-,i)|liy2, P(n) < C||en(-,t)|| LP(Q) < C, for almost every t G [0,T] (4.16) 

Thus, -B n satisfies (|4.3p ; that is 

ll-Dnlloo < C || Vf n || LOO (o i T;Ty 1 'P(n)) — Mlsn||oo < &■ (4-17) 

Furthermore, let v be the unique solution of (|4.ip corresponding to £. Then, by the linearity of Eq. 
(I4TTD and by (OB]l . 

IK(-,t) -v(-,t) \\ H 2(0) < C\£ n (-,t) -£(;t)\ 2 , for almost every t G [0,T]. 

Since £ n — > 77 strongly in L 2 (Q) in condition (|4.8p . it follows that 

||«n - v \\l 2 (0,T;H 2 (Q)) - C ll&» ~~ £\\l ~* °' 

which implies that 

B n = X Vv n -> xVv = £ strongly in L 2 (Q). (4.18) 

By (|4.14|) and (|4.17|) . the sequence {B n u n }^ =1 is bounded in L 2 (Q), whence there is a subsequence 
such that B n u n — >• 7r weakly in L 2 (Q) as n — > oo. By (|4.18p . there is a subsequence of {-B n K£Li 
such that B n — > B almost everywhere in Q. This together with (|4.15j) implies that B n u n — > Bu 
almost everywhere in Q. Therefore n = Bu, and 

B n u n — > Bu weakly in L (Q) asn^ oo. (4-19) 

Now, by (|4.12p . (|4.14p . and (|4.19p . we can pass to the limit as n — > oo in (|4.10p to get that u is a 
weak solution of (|4.2p corresponding to B and / in the sense of Definition [TJ Now, it only need to 
show that u G A(£) = ^(B). Actually, let U n = u n — u and F n = l u (f n — /). Then U n solves the 
system 

'dtU n = AU n -V- (B n U n ) - V • [{B n - B)u] + F n in Q, 
d v U n = on S, (4.20) 

u n {x,o) = o x g n. 
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Multiply the first equation of (|4.20p by U n , and integrate over £1 to give 

I- \U n \ 2 2 + \VU n \ 2 2 < C \\BnWlc \U n \ 2 2 + C ||«||^ |S n -B\l + C f F n Y n dx. 
at Jn 

It then follows from Gronwall's lemma that 

\U n (;T)\ 2 2 < e C ^ Bn ^ T (c ||«||^ ||B n - S||| + C / F„y„dxJ . (4.21) 

By (|4.12p . (|4.17p . and f|4.18[> . we see that the right-hand side of (|4.2ip tends to as n -> oo, and 
hence |C/ n (-,T)| 2 — > 0. Since u n (x,T) = for almost all 3; £ O, we get that u(x,T) = for almost 
all x G O. It then follows that u € ^(-6) = A(£). Therefore, A is upper semi-continuous. 

Finally, it remains to show that A(K) C K. Indeed, by the standard energy estimate we see 
that for any £ G K, each element u of A(£) satisfies 

Hull < p c ^+\\ b \\Io)( 1 + t ) (\ Vn \ _i_ in fll ^ 

which together with ([23]) in Proposition O and flHD leads to u G L°°(Q) n L°°(0, T; L p (fi)) and 

1 1 41 II -I- 111/ II <" <= ClK ° Ulr.1 

H U llL oo (0,T;LP(n)) + ll n lloo - e l U 0loo , 

where c\ is a positive constant independent of T, and kq is given by (|4.5p . If |«oloo — e _ClK0 which 
is exactly (II. 2p . then HuH^ < 1 and hence A(K) C -fC. Apply Kakutani's fixed point theorem to 
obtain at least one fixed point u of A; that is u G A(n). This it together with v = v u , the solution 
of (|4.ip with £ = u, gives the solution of (jl.ip . corresponding with some control / and u(x,T) = 0. 
This completes the proof. □ 
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